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Abstract 

In this paper, we use data collected from over 2000 non-residential electric vehicle supply 
equipments (EVSEs) located in Northern California for the year of 2013 to estimate the 
potential benefits of smart electric vehicle (EV) charging. We develop a smart charging 
framework to identify the benefits of non-residential EV charging to the load aggregators 
and the distribution grid. Using this extensive dataset, we aim to improve upon past studies 
focusing on the benefits of smart EV charging by relaxing the assumptions made in these 
studies regarding: (i) driving patterns, driver behavior and driver types; (ii) the scalability 
of a limited number of simulated vehicles to represent different load aggregation points in 
the power system with different customer characteristics; and (in) the charging profile of 
EVs. First, we study the benefits of EV aggregations behind-the-meter, where a time-of- 
use pricing schema is used to understand the benefits to the owner when EV aggregations 
shift load from high cost periods to lower cost periods. For the year of 2013, we show a 
reduction of up to 24.8% in the monthly bill is possible. Then, following a similar aggregation 
strategy, we show that EV aggregations decrease their contribution to the system peak load 
by approximately 40% when charging is controlled within arrival and departure times. Our 
results also show that it could be expected to shift approximately 0.25kWh (~2.8%) of 
energy per non-residential EV charging session from peak periods (12PM-6PM) to off-peak 
periods (after 6PM) in Northern California for the year of 2013. 


1 Introduction 

A recent analysis identifying the infrastructure and technology needs to meet California’s green¬ 
house gas (GHG) reduction goals for 2050 shows that the electrification of the transportation 
system plays a significant role in reaching these goals. In order to achieve the 80% reduction 
target in electrification, most of the direct fuel uses in buildings, transportation and industrial 
processes must be electrified. Among these, electrification of transportation yields the largest 
share of GHG reduction, where 70% of the vehicle miles traveled should be by electrically pow¬ 
ered vehicles [lj. A study by the Electric Power Research Institute (EPRI) [2] also suggests that 
electric vehicles will constitute a rather significant 35% of the total vehicles in the US by 2020. 
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This rapid growth in the electrification of transportation presents significant challenges as 
well as opportunities to the operation of today’s power system. When considered as inflexible 
loads, EVs will increase the current electricity demand significantly, intensifying the stress on 
the electric power system and pushing it closer to its limits. Research suggests that this is the 
case for low penetrations of EVs ESI El- However, when considered as flexible resources, where 
EV charging is controlled by direct or indirect strategies, EVs promote the reliable operation 
of the power grid 13 0 El, while also providing additional revenue streams that can be used 
towards the electrification of transportation [SEUS- This is particularly important considering 
the expected increase of renewable generation sources in the generation portfolio of many states 
in the U.S., as smart EV charging may provide the means to balance the intermittency of these 
resources. 

A number of recent studies aim to understand the adaptation needs of the existing operational 
control mechanisms to realize smart charging, and often propose novel planning and control ap¬ 
proaches. These approaches can be grouped into direct and indirect control approaches [?]. 
In direct control approaches, the control actions are realized without the vehicle owner in the 
control loop. Often, load aggregations are created to increase the size of the resource so it can 
offer economic benefits to the aggregator (8} [10]. In m, for example, the authors propose a 
direct load control strategy to provide vehicle-to-grid services for 3 different predefined mobility 
patterns. In [I2j, the authors conduct a simulation study for 3000 EVs parked at a municipal 
parking lot and evaluate the real-time performance of a direct control approach, which max¬ 
imizes the expected state of charge of the EV aggregation in the next time step subject to 
mobility constraints. In [13] . the authors develop an optimal direct control scheme based on 
global charging costs. The authors compare the proposed direct control scheme to the local 
scheduler in a simulation environment including 100-400 EVs. The arrival times of the EVs, the 
charging periods, and the initial energies of EVs are assumed to have a uniform distribution. 

The authors of [13 discuss various services that can be provided by electric vehicles, includ¬ 
ing peak shaving, regulation, voltage control, and reserves, and many studies have quantified 
the benefits of smart charging from various stakeholder perspectives [12113113. In [13, the 
authors demonstrate a proof of concept regulation case study. In G2, the authors estimate that 
smart charging will reduce the daily electricity costs of a plug-in hybrid EV by $0.23. They also 
identify daily profits for the individual driver when the charging of the vehicles can be regulated. 
The economic benefits of fleets that participate in specific markets have also been extensively 
studied. For example, in [T5] . 352 vehicles are used to estimate the economic potential of fleets 
when providing regulation up and down services using historical prices obtained from California 
Independent System Operator (ISO). In [IT], the authors use historical market data and charg¬ 
ing data collected from an EV located in a residential household to investigate financial savings 
and peak demand reduction. The authors conclude that the peak EV demand can be reduced 
by up to 56%. 

In this paper, we primarily focus on direct control approaches—in particular, centralized 
smart charging of EV aggregations—and we create a variety of case studies to investigate the 
potential benefits of smart charging to different stakeholders. To develop these case studies, we 
use data collected from over 2000 non-residential electric vehicle supply equipments (EVSEs) 
located throughout 190 zip code regions in Northern California spanning one year. To the best 
of our knowledge, this is the first study that uses such an extensive dataset on EV charging. 
First, we analyze over 580,000 charging sessions to investigate the trends in load flexibility and 
infrastructure use in the dataset. Next, we create virtual aggregation points (VAP) in which 
a combination of the EVSEs is assumed to be fed by the same distribution feeder. The VAPs 
mostly coincide with Pacific Gas and Electric Company’s (PG&E) sub-load aggregation points 
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(sub-LAPs). Additional details regarding this relationship is provided in Section [ 2 ] We in¬ 
troduce a smart charging framework to estimate the benefits of smart EV charging to various 
stakeholders in each VAP. As an initial case study, we investigate the potential benefits of EV 
aggregations operated under a single owner, where a time-of-use pricing scheme is used to es¬ 
timate economic benefits to the owner via shifting load from high cost periods to lower cost 
periods. Then, we create a case study where EV aggregations are used to decrease their current 
contribution to the system-level peak load. 

The motivation for this study is threefold: (*) Most of the work investigating the potential 
of smart charging of EVs is based on assumptions made regarding trip and customer character¬ 
istics. For example, in |18| . the authors use a fleet which includes commuter cars, family cars 
and taxis with predetermined departure and arrival locations randomly selected from a limited 
number of alternatives. In [19], the authors use data from driving surveys that reflect the driving 
behavior of people using internal combustion engine cars. They assume that the driving behav¬ 
ior of an EV owner will be similar to that of an internal combustion engine car owner. The 
dataset used in this study allows us to extract trip and customer characteristics, hence no such 
assumptions are needed on these characteristics, (ii) Often, a limited number of vehicles and 
mobility patterns are used in fleet-based studies to capture the most likely driving scenarios. For 
example, in CD, the authors develop a proof of concept strategy and show cost benefits for 50 
EVs with 3 different pre-defined mobility patterns. Although the exact number of EVs are not 
available in the dataset used in this study, the number of charging sessions (over 580,000) and 
the fact that these charging sessions are spread throughout the year ensure that a representative 
population of non-residential charging is studied, (iii) The individual charging profile of an EV 
is often represented by a typical constant-voltage, constant-current curve for certain battery 
chemistries, or more simply by a constant charging power [7J. For example, in [20] . the charging 
power is assumed to be fixed at 4.4kW, whereas in m, the authors use the charging profile of 
a typical lithium-ion battery pack obtained from [22]. The dataset used in this study includes 
time series of power measurements obtained every 15 minutes for each charging session. Hence, 
no assumptions are made on charging profiles of the vehicles, and individual charging data is 
available for each charging session. 

The remainder of the paper is organized as follows: Section [2] introduces the dataset and 
discusses the load flexibility and infrastructure use trends obtained from the dataset. Section [3] 
presents the smart charging strategy used in this study. Specifically, it discusses the framework 
and the underlying assumptions made when estimating the benefits to different stakeholders. 
Sections [4] and [5] describe the case studies completed in this research. Finally, Section [6] discusses 
the conclusions and opportunities for future work. 

2 Dataset 

The data used in this study is collected from individual EVSEs located in 16 different sub-LAPs 
in PG&E’s territory for the year of 2013. For each charging session (i.e. from plug-in to departure 
of an EV), the EVSEs report the start and end period of the charging, the plug-in and departure 
time stamps, the average power, and the maximum power (measured every 15 minutes), as well 
as the charging port type, the location (zip-code level), and the non-residential building category. 
Since the dataset includes the location information based on zip codes and some zip codes are 
fed by multiple sub-LAPs, we create virtual aggregation points (VAPs) for the zip codes that are 
fed by multiple sub-LAPs. This is done by combining the sub-LAPs’ identifiers. Table [l] presents 
the final list of VAPs in the dataset and total number of zip code regions forming each of these 
VAPs, the total number of charging sessions, and the average number of daily charging sessions 
in each VAP. Figure [l] depicts the centroids of the zip code regions forming the considered VAPs. 
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VAP 

Region 

# of zip 
code regions 

# of 
charging 
sessions 

# of 
charging 
sessions per 
day 

P2-SB 

Peninsula & South Bay 

7 

207501 

568.50 

SB 

South Bay 

21 

112250 

307.53 

SF 

San Francisco 

30 

72996 

199.99 

P2 

Peninsula 

17 

59252 

162.33 

EB 

East Bay 

27 

52700 

144.38 

EB-SB 

East Bay &; South Bay 

6 

16902 

46.31 

NB 

North Bay 

14 

12346 

33.82 

LP 

Los Padres 

8 

9035 

24.75 

CC 

Central Coast 

15 

8428 

23.09 

SA 

Sacramento Valley 

11 

7787 

21.33 

FG 

Geysers 

11 

7918 

21.69 

SA-SI 

Sacramento V. &; Sierra 

2 

7465 

20.45 

CC-P2 

Central Coast & Peninsula 

2 

6778 

18.57 

FG-NB 

Geysers & North Bay 

4 

3845 

10.53 

FI 

Fresno 

4 

377 

1.03 

NV 

North Valley 

1 

336 

0.92 

ST 

Stockton 

3 

244 

0.67 

FG-NC 

Geysers Sz North Coast 

1 

246 

0.67 

SI 

Sierra 

2 

181 

0.50 

SN 

San Joaquin 

1 

134 

0.37 

HB 

Humboldt 

1 

101 

0.28 

P2-SF 

Peninsula & San Francisco 

1 

73 

0.20 

NC 

North Coast 

1 

15 

0.04 


Table 1: VAPs used in this study 


In this study, we use data from VAPs with an average of 20 or more charging sessions per 
day. These VAPs are indicated in bold in Table [l] Figure [2] also shows the total number of 
charging sessions per month for each VAP used in this study. Over the course of 2013, the total 
number of charging sessions approximately doubles. 

Figure [3] shows the combined load profiles of VAPs for the second weeks of January and 
December. The impact of the growth in charging session is reflected on the daily load profile of 
the loads. Moreover, the peak non-residential EV load occurs between 9AM and 11AM, and it 
more than triples from January to December of 2013. 

2.1 Load Flexibility and Infrastructure Use 

To gain further insight into the dataset and to understand the distribution of charging sessions 
and the use of EVSEs in different regions, we analyze the charging sessions obtained from the 
VAPs marked in bold Table [T} The infrastructure use, I use , in each VAP is represented by the 
average number of charging sessions N sess i ons per EVSE and calculated for every business day 
of 2013. Formally: 
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Figure 1: Centroids of zip code regions forming the VAPs 
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where Nevse is the number of EVSEs. Figure [4] depicts the box plots of the infrastructure use 
within 2013 for all of the VAPs. For each month of 2013, a box plot is created to represent 
the distribution of the I use values calculated for every business day of the month. The median 
value of infrastructure use is marked with a red line in each box plot, and the boundaries of 
the box depict the 25th and 75th percentiles. The whiskers correspond to the 99th percentiles 
assuming the distributions per each month are normal. The median infrastructure use increases 
in all VAPs from 1.8 to 2.1 sessions per EVSE from January to December. We believe that this 
is due to the fact that the demand has increased faster than the number of EVSEs. 


In addition to the infrastructure use, we investigate the load flexibility in each VAP. The load 
flexibility depends on the charging duration d c harge and the overall duration of each charging 
session ^session- Formally, we define the load flexibility Zfl ex as the ratio of the duration that a 
car is plugged but not charging to the overall session duration: 


^flex — 


-d, 


charge 


( 2 ) 


Figure [5] depicts the load flexibility for all VAPs by month. As Figure [5] suggests, the load 
flexibility decreases slowly as the number of charging sessions per EVSE increases. Also, most 
of the distributions have a slight positive skew. The size of the box representing the 25th and 
75th percentiles is also decreasing with time, suggesting an increase in skewness. 


The load flexibility metric shows the charging duration relative to the session duration; how- 
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Figure 3: Load shapes for January and December for all the VAPs 
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Figure 5: The variation in load flexibility 


ever, it does not capture when the charging sessions occur. The start and end times of the 
charging sessions play a key role when estimating the benefits of EV aggregations to the power 
system. To put these two variables into perspective, we show a histogram of arrival (i.e. session 


start) and departure (i.e. session end) times in Figures 6a and 6b respectively. 


As can be seen in Figures 6a and 6b, most of the charging sessions start within the 7AM- 
10AM period and often end within the 5PM-7PM period. Considering these loads are currently 
uncontrolled (i.e. they immediately start charging when they are plugged in), they coincide with 
the typical working hours of a non-residential location. These figures suggest that employees or 
customers arrive in the morning and plug in their vehicles. Some leave around noon and come 
back, and most leave work between 4PM and 7PM. 


3 Smart Charging Strategy 

In this section, we introduce the proposed smart charging methodology. In particular, we de¬ 
scribe the general optimization strategy used to obtain the charging schedules for each charging 
session. 

The goal of the proposed smart charging framework is to reschedule the power time se¬ 
ries measured in discrete time slots [1,..., K] for any charging session in a population of EVs, 
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Figure 6: Arrival and departure time characteristics 

[P\, P 2 ,..., Pk\ such that an objective function is optimized. The objective function should 
capture the desired benefits from a stakeholder’s perspective. While rescheduling the charging, 
we would like to ensure that the order of the measurements in this time series is preserved. This 
is because the power that EVSEs draw is dependent on the state of charge (SOC) of the EV that 
is being charged, and keeping the order of the measurements accounts for this dependency. In 
addition, we assume that the charging is preemptive; that is, the charging tasks are interruptible 
without any decrease in the SOC of the EV. 

In a typical charging session, an EV starts charging when it is plugged in, and often the 
charging is complete before the vehicle departs. The smart charging framework proposed in this 
study is designed to move some of the charging to the slack time slots (i.e. the time slots where 
the vehicle is plugged in but the charging is completed). 

For the purposes of this paper, we discretize a day into 15-minute intervals. We define the 
time period for the optimization within a day as the time between the start time slot t s tart and 
the end time slot t en d . In this period, each charging session i has an arrival time slot denoted 
by to and a departure time slot t d . For each charging session, a column vector including the 
charging power time series can be created using the power measurements for every time slot in 
[ta\ t^']. If necessary, the time series is zero-padded to match the size of the optimization time 
period \t s tarti tend]- Hence, for each EV i, the power time series is given as follows: 

= (3) 

where K is the total number of time slots in [t s t ar t, tend]- Next, for each charging session i. we 
identify Qb) whose elements correspond to the j th non-zero element of Pd). The goal is 
to reschedule the time slots t!'^ in [ta\ t! d \ corresponding to without changing their order. 
We define Af-*) as the total number of non-zero power measurements in charging session i (i.e. 
total number of elements in Q w). 

To capture the precedence and the session duration constraints we proposed above, the 
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following formal constraints are introduced: 
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Vj G [1 ,Af«] 


( 4 ) 


The proposed constraints are constructed using a binary decision matrix to represent charging or 

(i) 

non-charging time slots within the optimization duration. In particular, for each element Q y - in 

QW, we create a binary vector x^V) that includes K binary decision variables. Each element in 

(i) 

this vector represents a candidate time slot at which Q y -' could be positioned. Hence, we define 

row vectors x^’-O \/i £ [l,jV] and Vj G The elements in these vectors are xl/"^ G {0,1} 

that are defined \/k G [ I, /V]. 


From these binary vectors x*-*V), we form a binary decision matrix for each charging 
session i G [1, iV]. In particular, the individual decision variables form the elements of the 
binary decision matrix X W as follows: 


JfW = 


t (M) 

x 1 

T (01) 

x K 

(i,M W) 

1 

(i,Af ( 0) 

• ’ X K 


( 5 ) 


Finally, we write the variables in the constraints given in Q using the binary decision variable 
as follows: 


fW = X (i) 0, where O 


1 

2 


it' 


( 6 ) 


The aggregate power vector for the YAP AP^ = YliL 0 (pV)j 


for the day d is given as follows: 




T 

rx (1) i 

A p( d ) = 



X( 2 ) 


-Q {N) . 


l( W) . 


For each case study proposed in this paper, we build on the general optimization frame¬ 
work described above, identify the objective functions to capture the benefits from each stake¬ 
holder’s perspective and introduce additional constraints when necessary. We use the Gurobi 
optimizer |23| to solve the optimization problems formulated for each case study. Due to the size 
of the optimization problem for certain VAPs and the number of times the optimization problem 
is solved to obtain values to estimate benefits for the year of 2013, a proved optimal solution is 
expected to be hard to reach within a reasonable time frame. For these reasons, we alter the 
optimality criteria by controlling the relative gap between a feasible integer solution and the 
general optimal solution. We set this optimality criteria to 5% and allow early termination once 
a feasible solution is found. 
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4 Charging Infrastructure Owner’s Perspective 


In the first case study, our goal is to capture and maximize the benefits of smart charging from 
an EV charging service provider’s perspective. Currently, each charging meter is independently 
owned by the building owner, and the consumption is billed to the building owner as part of 
the building’s monthly bill. However, in our work, we focus only on the load resulting from EV 
charging, i.e. decoupled from other loads, but aggregated over VAPs formed based on sub-LAPs. 
This corresponds to the situation in which the charging stations within each VAP are combined 
and operated under a single owner or an aggregator and the owner is charged according to a time 
of use (TOU) tariff structure, where shifting load from high cost periods to lower cost periods 
can offer some benefits to the owner. Although the current VAPs are created based on sub- 
LAPs, the current scale of the charging infrastructure and the number of charging sessions can 
easily represent a large parking structure or a campus in the future, where the EV aggregation 
is behind a single meter and the non-EV load is relatively steady. 


4.1 Problem Formulation 

In a typical TOU rate structure, there are two separate charges forming the monthly bill: the 
energy charges and the demand charges. The energy charges are calculated based on the amount 
of energy consumed over given time periods of the day using the corresponding hourly TOU en¬ 
ergy rate. The demand charges are calculated based on the maximum power demand for specific 
time periods of the day over the course of the billing period. At the end of each billing period, 
the maximum demand values for the specified periods are multiplied by the demand charge rates 
and added to the overall energy charge. 

In order to model a similar rate structure in the proposed smart charging framework, we 
define EC ^ as the energy charge for day d of a month with D days (i.e. d 6 [1,..., D ]). Then, 
we define DCh as the demand charges for each time period h of the day of any month. For 
example, in PG&E’s E-19 TOU rate structure, for winter billing periods, the demand charges 
are calculated based on 2 time periods paid-peak (i.e. 8:30AM-12:00PM k 6:00PM-09:30PM) 
and off-peak (i.e. 09:30PM-08:30AM) [23] ■ Formally, the monthly bill for the owner is therefore 
given by: 


f(DC h , EC «) =J2 DC h + J2 ECid) (8) 

V/i Vd 

The energy charges EC ^ can easily be incorporated into the proposed daily optimization 
routine. Defining ER as a column vector reflecting the price of energy for each time slot j, 
ECW for any day d in a billing period is given by: 

EC (d) = AP^ER (9) 

For time period h within day d, a subset of the entire daily aggregate power vector AP ^ is 
needed and is referred to as AP^ d \ 

In order to minimize the cost function given in (|8]), the maximum demand for the daily 
time periods h must be accurately known beforehand for the entire month. However, in a real 
life scenario, this is not a valid assumption. To incorporate demand charges into the proposed 
daily smart charging framework, we therefore propose the following strategy for the owner: for 
each day d , we define the peak aggregate power values for each period h as AP^ ak h . Since the 
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historic AP pea ^values for each day in [1,..., d — 1] are available to the main scheduler, we can 
define the maximum of the historic AP pea k,h values until d — 1 as follows: 

= rnax(AP^ t h ,.... AP^k) (10) 

Using the above definition, the monthly demand charges can be calculated at the end of the 
month based on AP^ h and the demand rates DRh for each period as: 

OO h = AP^ h DR h ( 11 ) 

As we move from one day to the next, we try to limit the demand charges based on the 
maximum daily demands occurred up to the current day. At the beginning of the billing period, 
we start with no knowledge of the historical peak values, and we keep track of the maximum 
historical value up to day d. This strategy can be represented by incorporating the maximum 
value of the peak values AP^ ax , for time period h and day d as decision variables into the 
following optimization problem: 


minimize 

X(i),AP W , 

’ max,ri 


AP 


(d) 


max.h 


DR h + EC^ d) 


subject to Q and the following additional constraints: 


AP 


(d~ 1) 


max,h — 

M 


< AP, 


(d) 


max,h 




A p(“) < A p 
Ar h — Ar max,h 


Vh E [1, TP] 


( 12 ) 


Note that with (12), we ensure that the current maximum AP^ ax . is more than or equal to 

the maximum historical value AP for period h. By definition, this allows for the tracking 
of the maximum value up to that day. In addition, these maximum values set the day based 
on which the demand charges will be calculated. If none of the current peak values exceeds the 
historical maximum values, the demand charges for each period h are not set by the current day d. 


4.2 Case Study 

For the purposes of this paper, we use the demand and energy rates from PG$E’s E-19 TOU rate 
structure [24]. The E-19 rate structure gives the owner the option to manage their electric costs 
by shifting load from high cost periods to lower cost periods. Detailed information on E-19 is 
given in Table [2] The summer period starts with May 1st and ends October 31st, and the winter 
period includes the remaining months of the year. This rate is for non-residential customers in 
PG&E’s territory with highest demand exceeding 499 kW for three consecutive months. 

To evaluate the benefits of smart charging when the EV aggregation has a single bill cal¬ 
culated on a TOU tariff, we first calculate the current bill under this tariff but without smart 
charging. Then, we use the proposed optimization strategy to schedule the loads in a way that 
minimizes the customer’s monthly bills, and we report each monthly bill calculated for each 
VAP and the contributions from energy and demand charges in the bill. 

4.3 Results 

Figure [7] shows the sum of monthly bills calculated in dollars for all of the VAPs. For each month, 
the left bar shows the current bill, and the right bar shows the optimized bill for the month. 
It is obvious that the difference between the summer and winter rates impacts the aggregate 
monthly bill. The increase within the winter and the summer period is due to the increase in 
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Demand Charges $/kW Time Period 


Max. Peak Demand Summer 

Max. Part-Peak Demand Summer 

Max. Demand Summer 

Max. Part-Peak Demand Winter 
Max. Demand Winter 

Energy Charges 

$19.71253 

$4.07 

$12.56 

$0.21 

$12.56 

$/kWh 

12:00PM-6:00PM 

8:30AM-12:00PM 

& 

6:00PM-09:30PM 
Any time 
8:30AM-09:30PM 
Any time 

Time Period 

Peak Summer 

$0.16253 

12:00PM-6:00PM 

8:30AM-12:00PM 

Part-Peak Summer 

$0.11156 

& 

6:00PM-09:30PM 

Off-Peak Summer 

$0.07818 

09:30PM-08:30AM 

Part-Peak Winter 

$0.10479 

08:30AM-09:30PM 

Off-Peak Winter 

$0.08200 

09:30PM-08:30AM 


Table 2: E-19 rate structure [2]ff 


the number of charging sessions over the year. 


Figures 8a and |8b| show the total energy and demand charges, respectively, over all LAPs. 
The cumulative energy charges increase slightly for the summer months when using smart charg¬ 
ing, whereas there is a significant drop in the demand charges. This suggests that the peak load 
of the EVs is shifted from the morning partial-peak period to the peak-period. This shift is still 
beneficial because the increase in the energy charges is significantly lower than the decrease in 
the demand charges. 


The cumulative load shapes given in Figure [3] and the arrival and departure time histograms 
given in Figures [6a] and [6b] support these results. These figures suggest that energy charges 
increase because a large portion of the EV charging sessions end (i.e. the charger is unplugged) 
before the system peak period ends. Thus, when coupled with the higher part-peak demand 
rates, the optimization converges to a result in which the load is shifted from the EV load peak 
period (9AM-11AM) to the system peak period (12PM-6PM). 

The results given in Table [3] provide further insight into the results depicted in Figures [7] 


for winter and summer months. Then, we report on average bill reduction per session during 
these periods. The values range between 0.13 and 0.65 dollars among all VAPs. Overall, we 
find that the rate structure in the summer periods yields to more reductions per session than 
the rates in winter months, with the exception of the Sacramento Valley (SA) VAP. We also 
report on the total percent bill reduction and we break down this percentage into contributions 
from demand charges and energy charges. We observe that average percent bill reductions range 
between 8.63% and 24.77%. Even though the average reduction per session values are mostly 
higher during summer months, the relative cost reduction in monthly bills for individual VAPs 
varies less. This is due to high overall costs in the summer months. 


8a and |8b| Specifically, we reflect on the average monthly bill before and after optimization 


Figures [9a] and [9b] depict the relationship between the reduction in demand charges and the 
number of charging sessions in each VAP per month. Specifically, in Figure [9a] we examine the 
decrease in demand charges in dollars. We observe a linear trend: as the number of sessions 
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Figure 7: Monthly bills calculated with E-19. The left bar for each month shows the current 
bill, and the right bar shows the optimized bill. 


VAP 

Period 

Bill [dollars] 

Reduc¬ 

tion 

Reduction [%] 



Cur¬ 

rent 

Opti¬ 

mized 

[dollars 

/session] 

DC 

EC 

Total 

P2-SB 

Summer 

63001 

50395 

0.65 

20.86% 

-0.85% 

20.01% 

Winter 

29603 

22575 

0.46 

23.41% 

0.33% 

23.74% 

EB-SB 

Summer 

4588 

3788 

0.52 

16.96% 

0.49% 

17.45% 

Winter 

2092 

1724 

0.28 

17.23% 

0.36% 

17.59% 

SA-SI 

Summer 

1645 

1413 

0.36 

13.80% 

0.30% 

14.10% 

Winter 

828 

752 

0.13 

9.06% 

0.12% 

9.18% 

CC 

Summer 

2365 

2178 

0.24 

7.34% 

0.57% 

7.91% 

Winter 

1037 

896 

0.22 

13.31% 

0.29% 

13.60% 

EB 

Summer 

12033 

10003 

0.41 

16.44% 

0.43% 

16.87% 

Winter 

5874 

4868 

0.26 

16.66% 

0.47% 

17.13% 

FG 

Summer 

1803 

1568 

0.33 

11.98% 

1.05% 

13.03% 

Winter 

920 

807 

0.18 

11.82% 

0.46% 

12.28% 

LP 

Summer 

2370 

2135 

0.29 

9.37% 

0.55% 

9.92% 

Winter 

1141 

1002 

0.20 

11.88% 

0.30% 

12.18% 

NB 

Summer 

3136 

2865 

0.23 

8.16% 

0.49% 

8.64% 

Winter 

1391 

1271 

0.13 

8.48% 

0.22% 

8.63% 

P2 

Summer 

16795 

14171 

0.48 

16.13% 

-0.51% 

15.62% 

Winter 

8567 

7010 

0.34 

17.98% 

0.20% 

18.17% 

SA 

Summer 

2313 

1991 

0.45 

13.88% 

0.04% 

13.92% 

Winter 

1215 

914 

0.52 

24.76% 

0.01% 

24.77% 

SB 

Summer 

32911 

27439 

0.53 

17.72% 

-1.09% 

16.63% 

Winter 

15645 

12602 

0.37 

19.34% 

0.11% 

19.45% 

SF 

Summer 

17679 

14224 

0.51 

18.07% 

1.47% 

19.54% 

Winter 

8591 

7046 

0.28 

17.10% 

0.88% 

17.98% 


Table 3: Average results based on summer and winter month rates in E-19 
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Demand Charges Energy Charges 

[Thousands of Dollars] [Thousands of Dollars] 



Jan Feb Mar Apr May June July Aug Sep Oct Nov Dec 

Month of 2013 

(a) Monthly energy charges calculated with E-19 
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(b) Monthly demand charges calculated with E-19 

Figure 8: Decomposition of monthly bills to energy and demand charges. In each figure, the left 
bar shows the current charges, and the right bar shows the optimized charges for each month. 
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Figure 9: Demand Charge Reduction by Session Size 
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per month rises, the reduction in demand charges increases linearly. Given the current load 
flexibility and arrival and departure times, this is expected because most of the EVs contribute 
to the peak load of the EV aggregation. In Figure 9b, we look at the percent reduction in 
demand charges. For up to 2000 charging sessions per month (indicated by a red dashed line in 
Figure [9b] ), there is no clear separation between the winter and summer months and, for a given 
number of sessions, the demand charge reduction values vary. Beyond this point, we can see 
a clear separation between the winter and summer months, and the demand charge reduction 
values show less variance. 


The relative decrease in the summer months is less than the relative decrease in the winter 
months. We believe that this is due to the time of the peak EV load, the arrival and depar¬ 
ture patterns of the EVs and the corresponding rate structure. In particular, the peak EV load 
coincides with the part-peak rate period, and most of the EVs depart before the system peak 
period (12PM-6PM) is over. The system peak period has a separate and higher demand rate in 
the summer months (detailed in Table [2J. This limits the smart charging framework’s ability to 
move the EV loads from part-peak period to system peak period. The winter rates we use in this 
study do not include a separate demand rate for the system peak period; rather, the part-peak 
period extends from 8:30AM-09:30PM. This makes it possible to manage the EV peak load in 
a more effective way. 


5 Distribution System Operator’s Perspective 

In the second case study, we evaluate the potential benefits that smart charging can offer to 
distribution system operators (DSOs). The motivation behind this case study is to investigate 
the potential of each charging session to decrease its contribution to the peak system demand via 
smart charging. We first quantify the percentage of peak load shed during the system peak load 
period (12AM-6PM). We then quantify the amount of energy that is shifted outside the peak 
period by the EV load aggregation for each month of 2013. Finally, we report on the amount of 
energy that can be expected to be moved outside of the system peak period per charging session. 


5.1 Problem Formulation 


To realize peak shaving, we propose to develop a two-stage optimization. In the first stage, we 
minimize a bound on the aggregate power consumed by the EVSEs within a VAP during the 
pre-defined peak period (12AM-6PM) only. We simplify refer to the pre-defined peak period as 
pp, and to simplify the notation introduced earlier, we refer to the aggregate power vector within 
the peak period as APpp . To implement the initial stage optimization, we define APf^ nd as 

a decision variable to represent the proposed bound on the APpp . Then, in the second stage, 
using the optimal bound as a constraint, we minimize the total energy consumed within the peak 
period. This implicitly ensures that the energy bill for the customer is decreased or unchanged 
based on a typical TOU tariff. The first part of the optimization can be written as: 


minimize 


AP, 


(d) 


X(i) AP, (d) , 

’ bound,,pp 


bound,pp 

subject to Q and the following additional constraints: 


a p ( d ) < a 
^ r w — bound,pp 


(13) 
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Figure 10: Distribution of percent peak shed for all the VAPs 
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Then, using the optimal AP^ o ^ ndpp values obtained in the first stage , we can form 

the second stage as follows: 


minimize 

x(>) 


£ APf 

VfcC pp 


subject to Q and the following additional constraints: 


A p( d ) < A p( d ) 
- f±r pp — 1 ‘ bound,pp 


(14) 


5.2 Case Study 

The motivation behind our second case study is to evaluate the potential of EV aggregations 
to decrease their contribution to the system peak load via smart charging. As the arrival and 
departure time histograms given in Figures [6a] and [6b| suggest, the amount of energy that can 
be moved outside of the peak period is expected to be low, mostly because most non-residential 
EV sessions end before the system peak period is over. However, there is potential in using 
smart charging and exploiting the inherent flexibility in each charging session to decrease the 
contribution of EVs to the system peak load. 

To demonstrate and quantify this potential, we calculated optimal schedules for each VAP- 
level aggregation using the optimization strategy described in the above section, and obtained 
percentage of peak shed values and the total energy moved outside of the peak period for every 
day in each month of 2013. 


5.3 Results 


Figure 10 shows the box plots created using daily peak shed values for each month of 2013. The 
percentage of peak shed for each day d is defined as: 


/p( d ) 

%peakshed ^ =- bound,pp (15) 

max(APpp ) 

The red lines denote the median value of the distribution, the box boundaries are the 25th 
and 75th percentiles and the whiskers denote the 1st and 99th percentiles, assuming the distri¬ 
butions per each month are normal. The outliers outside the whiskers’ boundaries are marked 
with points. As expected, the smart charging significantly reduces the peak EV load during the 
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system peak period. The median values for all of the months range between 30 and 42%. A 
decrease in the peak shaving potential and a slight decrease in the variation of the distributions 
over the course of 12 months are also apparent in Figure 10 This can be explained by the 


increase in the number of charging sessions per EVSE and the related decrease in the variation 
of available flexibility. 


Figure 11 depicts the distribution of the average energy moved outside of the peak period 


per charging session for all of the VAPs estimated every day of the month. The median value 
over 2013 is approximately 0.25kWh per charging session, which corresponds to ~2.8% of the 
average energy put during each charging session. 


6 Conclusions and Future Work 

In this paper, we quantify the potential benefits of smart charging to different stakeholders using 
data collected from over 2000 non-residential electric vehicle supply equipment (EVSEs) located 
throughout 190 zip code regions in Northern California. We created virtual aggregation points 
(VAP) in which the aggregate power consumption of a selected population of EVSEs is assumed 
to be managed via individual charging control at each EVSE. We developed and used a smart 
charging framework to estimate the benefits of EV smart charging to different stakeholders: a 
single owner/an aggregator of behind-the-meter EVSEs (i.e. aggregators) and distribution sys¬ 
tem operators. 

In our first case study, we investigated the potential benefits of behind-the-meter EV aggre¬ 
gations. The aggregate load is re-scheduled using a TOU rate structure. Our results suggest 
that up to 24.8% decrease in the aggregate monthly bill per VAP is possible. In all months, this 
reduction is due to a corresponding decrease in demand charges in the monthly bill: we observed 
that decreases in energy charges are contributing by up to 1.5% to the overall decrease, whereas 
the demand charges contribute up to 24.7%. 

In our second case study, we used the EV aggregations to decrease their contribution to the 
system-level peak load. We have observed median peak shed values around 30%-42% for each 
month. In addition, we have quantified the amount of energy that can be shifted outside the 
peak period per charging session over the course of 2013, and found the median value to be 
approximately 0.25kWh/session (~2.8% of the average energy put in each session). 
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In the future, we would like to investigate the impact of different non-residential customer 
categories (e.g., retail vs. workplace) within each VAP to similar metrics calculated in this study 
and identify suitable grid services for these customer categories. In addition, we would like to 
expand the current smart charging framework and develop control algorithms for workplace 
charging that use variable charging rates. We also would like to study the impacts of smart 
non-residential EV charging to the overall system load, in particular when the system level solar 
generation is expected to cause over-generation and ramping problems in the grid. 
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